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The  Euler-Bernoulli  Beam  Equation 
with  Boundary  Energy  Dissipation 


G.  Chen,  S.G.  Krantz*,  O.U.  Ma ,  C.E.  Wayne 
Department  of  Mathematics 
The  Pennsylvania  State  University 
University  Park,  Pennsylvania 

H.H.  West 

Department  of  Civil  Engineering 
Pennsylvania  State  University 
University  Park,  Pennsylvania 

INTRODUCTION 


Many  problems  In  structural  dynamics  Involve  stabilizing  the  elastic 
energy  of  partial  differential  equations  such  as  the  Euler-Bernoulli 
beam  equation  by  boundary  conditions.  Exponential  stability  is  a  very 
desirable  property  for  such  elastic  systems.  The  energy  multiplier 
method  f-t+r— {-2j-r-fTj‘"has  been  successfully  applied  by  several  people  to 
establish  exponential  stability  for  various  PDEs  and  boundary  conditions. 
However,  it  has  also  been  found ^£2]  that  for  certain  boundary  conditions 
the  energy  multiplier  method  Is  not  effective  in  proving  the  exponential 
stability  property. 

A  recent  theorem  of  F.L.  Huang  [4]  introduces  a  frequency  domain 
method  to  study  such  exponential  decay  problems.  In  this  paper,  we 
derive  estimates  of  the  resolvent  operator  on  the  imaginary  axis  and 
apply  Huang's  theorem  to  establish  an  exponential  decay  result  for  an 
Euler-Bernoulli  beam  with  rate  control  of  the  bending  moment  only.  We 
also  derive  asymptotic  limits  of  elgenfrequencies .  which  was  also  done 
earlier  by  P.  Rideau.[8]?  Finally,  we  indicate  the  realizability  of 
these  boundary  feedback  stabilization  schemes  by  illustrating  some 
mechanical  designs  of  passive  damping  devices. 
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§1.  BACKGROUND  AND  MOTIVATION 

In  this  paper,  we  consider  the  following  uniform  Euler-Bernoul 1 i  beam 
equation  with  dissipative  boundary  conditions 


«ytt(x.t)>Elyxxxx<x.t)  -  0. 

y(0, t)  *  0, 
yx(0,t)  *  0. 

~£Iyxxx(1-t)  ■  -k?ytd.t)  . 
-Elyxx(l,t)  *  k|yxt(l,t)  . 
(y(*,o),yt(*,Q)l  -  (y0(x),y1(x)), 


0  <  x  <  1 , 

kj  e  R, 

kj  e  i, 

0  <  x  <  1 . 


(1.1) 


where  ■  denotes  the  mass  density  per  unit  length,  El  is  the  flexural 
rigidity  coefficient,  and  the  following  variables  have  engineering 
meanings : 


y  =  vertical  displacement,  yt  =  velocity 
yx  *  rotation,  yxt  -  angular  velocity 
-Elyxx  -  bending  moment 
-EI*xxx  *  shear 

at  a  point  x,  at  time  t. 

From  now  on.  when  we  write  equation  (l.l.j).  for  example,  we  mean 
the  Jth  equation  in  (1.1). 

The  above  equation  and  conditions  are  intended  to  serve  as  a  simple 
mathematical  model  for  the  mast  control  system  In  NASA's  COFS  (Control 

of  Flexible  Structures)  Program.  See  Figure  1.  A  long  flexible  mast  60 
meters  in  length  is  clamped  at  its  base  on  a  space  shuttle.  The  mast  is 
formed  with  54  bays  but  can  be  idealized  as  a  continuous  uniform  beam. 

At  the  very  end  of  the  mast,  a  CMG  (control  moment  gyro)  is  placed  which 

can  apply  bending  and  torsion  rate  control  to  the  mast  according  to 
sensor  feedback. 

Boundary  conditions  (1.1.2)  and  (1.1.3)  signify  that  the  beam  is 
clamped,  at  the  left  end,  x  =  0,  while  boundary  conditions  (1.1.4)  and 
(1.1.5)  at  the  right  end,  x  *  1,  respectively,  signify 

'shear  ( -El yxjtx )  is  proportional  to  velocity  ( y t ) 

bending  moment  <-EIyxx)  is  negatively  proportional  to  angular 
.velocity  (yxt) 

Thus  the  rate  feedback  laws  (1.1.4)  and  (1.1.5)  reflect  some  basic 
features  of  the  CMG  mast  control  system  in  COFS. 
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mast 


space 


shuttle 


j  Figure  1  Spacecraft  mast  control  experiment 

The  elastic  energy  of  vibration,  E { t ) ,  at  time  t,  for  system 
(  (1.1)  is  given  by 

! 

Eft)  «  \  /£[«y*(x.t)  *  EIy*x(x.t)]dx. 

Note  that  In  (1.1)  we  have  already  normalized  the  beam  length  to  1. 

The  qualitative  behavior  of  (1.1)  has  been  studied  in  an  earlier 

paper  [2).  There  It  Is  shown  that  if  k*  >  0.  k|  >  0  in  (1.1.4)  and 

(1.1.3),  respectively,  then  the  energy  of  vibration  of  the  beam  decays 
uniformly  exponentially: 


E(t)  <  Ke*"tE(0)  (1.2) 

for  some  K.m  >  0  uniformly  for  all  initial  conditions  (yQ(x) . y } (x ) ) . 

Therefore  the  flexible  mast  systen  can  be  controlled  and  stabilized. 

The  proof  of  the  above  in  [2]  was  accomplished  by  the  use  of  energy 
multipliers  and  the  construction  of  a  Liapounov  functional. 

Nevertheless,  a  major  mathematical  question  remained  unresolved  in  [2): 


[Q]  "Does  the  uniform  exponential  decay  property  (1.2)  hold  under  the  Accession  For 

NTIS  GFA&I 


2  2 

assumption  of  kj  *0,  k£  >  0?" 


DTI C  TAR  TJ 

Unannounced  Q 


This  question  is  of  considerable  mathematical  interest  because  the  |  Justification- 
feedback  scheme  using  bending  moment  only  is  simple  and  attractive.  | _ _ — 
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For  a  long  tine,  we  have  conjectured  that  the  answer  to  [Q]  Is 
affirmative,  as  asymptotic  elgenfrequency  estimates  obtained  In  [6]  (and 
§3)  have  so  suggested:  Let  A  denote  the  Infinitesimal  generator  of 

2 

the  C0-senigroup  corresponding  to  (1.1)  with  kj  *  0  and  k|  >  0.  and 
let  o(A)  denote  the  spectrum  of  A.  Then  there  exists  0  >  0  such  that 

Re  X  <  -0  <  0  for  all  X  e  o(A).  (1.3) 

Nevertheless,  it  Is  well  known  [4]  that  the  following  "theorem"  is 
false. 

"Let  A  generate  a  Cp-semigroup  and 

sup(Re  X|X  €  0(A))  <  -0  <  0  (1.4) 

for  sone  0  >  0.  Then  the  Cg-semigroup  Is  exponentially  stable: 

||exp(tA)||  <  Mc_Ait  for  some  M  >  1,  >  0".  (1.5) 

Therefore,  knowing  (1.4)  alone  is  not  sufficient  to  confirm  (1.5).  This 
statement  remains  false  even  if  we  assume  additionally  that  A  has  a 
compact  resolvent. 

We  have  repeatedly  tried  to  refine  the  energy  multiplier  technique 
used  in  [2]  to  establish  (1.5)  without  much  success,  no  matter  how  many 
different  and  elaborate  multipliers  were  constructed.  There  always  are 
boundary  terms  which  cannot  be  absorbed  by  terms  in  the  dissipative 
boundary  condition. 

A  recent  theorem  by  F.L.  Huang  offers  an  Important  direct  method 
for  proving  exponential  stability: 


THEOREM  1  (F.L.  Huang  f 4 1 


Let  exp(tA)  be  a  Cn-semigroup  in  a  Hilbert  space  satisfying 


l| exp ( t A ) ||  <  Bg,  t  >  0.  for  some  BQ  >  0.  (1.6) 


Then  exp(tA)  is  exponentially  stable  if  and  only  if 


(lu|u  €  R)  C  p( A),  the  resolvent  set  of  A;  and 


(1.7) 
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are  satisfied.  Q 

Huang's  r.heoren  effects  a  frequency  domain  method  to  proving 

exponential  decay  properties.  As  mentioned  earlier,  the  energy 
Multiplier  Method,  which  corresponds  to  a  time  domain  method,  has  not 

2  2 

been  successful  for  the  case  kj  *  0,  k:;  >  o. 

Therefore  the  work  is  to  obtain  bounds  on  the  resolvent  operator 

(lw-A)*1.  Here  we  accoeplish  this  by  carrying  out  a  careful  analysis 
on  the  eigenfunctions  and  eigenfrequencies  of  the  operator  A.  This  is 
done  in  §2. 

Associated  with  [Q]  is  the  question  of  the  asymptotic  distribution 
pattern  of  eigenfrequencies,  as  nuaerical  study  in  (2)  suggests  that  a 
"structural  damping"  phenomenon  is  present  at  low  frequencies.  Does  it 
also  appear  at  high  frequencies?  This  is  answered  in  §3.  (We  must  state 
that  the  work  and  nuaerical  verification  was  done  ahead  of  us  by  P. 

Rldeau  in  his  recent  thesis  [8]). 

In  §4.  we  present  Mechanical  designs  of  devices  satisfying  damping 
boundary  conditions  (1.1.3)  and  (1.1.4)  to  indicate  the  realizability 
of  the  feedback  stabilization  scheme  using  passive  dampers. 

2 

Notations :  We  use  ||  ||  to  denote  the  t  (0.1)  norm.  We  define  the 

Sobolev  space 

Hk  •  Hk(0.1)  -  (f :  [0.1J-*R)||f|j2  «  2  /o)f  (J)(x)|2dx  <  «>} ,  k  6  N. 

Hk(0 . 1 ) 

Also,  we  let 

Hq  ■  Hq ( 0 . 1 )  ■  (f|  f  €  H2(0.1),  f(0)  ■  f’  (0)  -  0). 

The  underlying  Hilbert  space  H  for  the  PDE  (1.1)  is 

»  -  Hq(0 , 1 )  X  r2(0.1)  ((f.g)lll  (f.g)ll^  •  /g[EI  |  f"  ( x )  |  2  *m|  g(  x)  |  2  )dx<<» } . 

whose  norm  square  is  the  elastic  energy. 

The  unbounded  linear  operator  A  associated  with  (1.1)  is  given  by 


with  doaaln 


o(A)  -  <<r.g)  €  h4*h2|  -Eif"( i )“-kjg( i ) . -Eif  ( i )*k|g‘  ( i ) , f (0) >ra  (o)  -o 


§2 .  ESTIMATION  OK  THU  RESOLVENT  OPERATOR  ON  THE  IMAGINARY  AXIS . 
EXPONENT  I AI.  DECAY  OF  SOLUTIONS. 

Consider  the  resolvent  equation:  Given  (f.g)  6  H  and  X  6  C,  find 
(«x,vx)  6  0(A)  such  that 


This  amounts  to  solving  the  following  boundary  value  problem  for  wx: 

a4w£4)(x)  ♦  X2wx(x)  -  - [xf (x ) «g(x)  | ,  x  e  (0.1) 
wx(0)  -  o 

wx  (0)  -  0 

*«J(i)  -  xk2wx(i)  -  itjf(i) 

w"(l)  ♦  Xk2w^(l)  -  -k|f  (1) 

where 

k2  «  kf ( El  )_1 .  k|  ■  k| ( El )  1 .  (2.3) 

Once  wx  is  found  we  obtain 

vx(x)  -  Xwx(x)  *  f(x)  (2.4) 

To  simplify  notation,  from  now  on.  unless  otherwise  specifically 
mentioned,  we  set  a4  -  1  in  (2.2.1)  and  write  kj ,  k2  for  kj .  k2 , 
respectively. 

The  main  work  in  this  section  is  to  prove  estimate  (1.8).  l.e..  to 
show  the  existence  of  some  Bj  >  0  such  that 


Eultr-Bernoulli  Beam  Equation 

/on“x(x)|2*|vx(x)|2]dx  <  B,/3[|f"(x)|2+|g(x)|2]dx  (2.5) 

for  all  X  »  lu,  u  e  K  and  all  (f,g)  6  K. 


LEMMA  2  A-1  exists  and  is  a  compact  operator  on  H.  Furthermore, 

o(A)  consists  entirely  of  Isolated  eigenvalues. 


Proof :  Let  X  -  0  in  (2.2),  We  see  that  wQ  in  (2.2)  is  obtained  by 
Integrating  four  times: 


"o(x>  -  r0f;f:  f0*  *<*!><«, <*2*3^ 

If'  Jo28«i>d<id^2  -  kif,<3>] 


x2 

T 


'  (£  -  rt[Jo  B(t,d< 


and 


vQ(x)  -  f (x ) 

Thus  A'1  exists  and  maps  M  into  K4 (0,1)  *  Hq (0,1).  Therefore  A'3 

is  compact.  The  rest  of  the  lemma  follows  from  Theorem  6.29  in  (6, 
Chapter  3) .  0 


LEMMA  3  The  resolvent  estimate  (2.5)  holds  for  X  »  iu,  u  6  R, 
provided  that  |X|  is  sufficiently  large. 

Proof :  For  simplicity,  let  us  write  (w,v)  for  (w^.v^)  when  no 

ambiguities  will  occur. 

Let 
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X 


i«  ■  Ijj2  , 


17  *  0. 
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n  in  -v  -in 

i  •  (»J3-In2kf )en  -(i7J3*iiJ2kf)e1’J  -('l3*i*i2k2)e~’'  ( in3- lnkf ) e- 1 17 

[n2*lr,3kl)ev  (-t,2Vk|)e1’7  (^2-lTj3k|  )e_T»  (-v2*v3  k|)e-1’' 


If  N*1  exists  for  v  sufficiently  large,  then 


(2.12) 


(2.13) 


»1  o 

He  now  begin  the  estinatlon  of  |  w*‘ { x )  I  dx.  The  work  below  may 

"  0 

••(■tedious,  but  the  idea  is  rather  simple.  The  main  observation  is 
th»t  the  dominant  terns  in  wp(x)  and  w(x)  do  not  satisfy  the  bounds 

f1  |w"(x)|2dx  <  cf1  (|f"(x)i2  *  |g(x)|2)dx 
•*0  p  J0 

f1  | w" (x) | 2dx  i  cf1  [ | f" (x) | 2  *  | g( x ) | 2 ) dx 
•  n  •  n 


f*r  |X|  large.  However,  in  (2.9),  those  dominant  terns  cancel . 
lMvlng  w(x)  with  smaller  terms  which  are  bounded  by  <J(||f"||  ♦  ||gj|) 


l»t  Step  Estinatlon  of  W'(x). 


ft(a  (2.8). 


■J(x)  *  -  1J*1C»inh  >;(x-f)*sin  ij(x- (  ) )  ( l*>2f  ( (  )  *g(  (  )  ]df  (2.14) 
X 

'-if  «~l(sinh  ij(x-{  )+sln  tj(x  (  )  ]g(f  )d(  (integration  by  parts) 

*•  n 


if  ij_1(sinh  n(x-f)-sln  ij(x-f )]if"(f )df 

ft 


'  •:  -V*'  '•V 
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-  -  I  -f1/*  slnh  ►  Olv^llir'lhllell) 

«  -  \  ’j'1J^en(x't)[if"(f  )^K(f)]^  ♦  0(T|-1  [||  f"|Mlg||]) 

-  -  i  ij'1  e,,xJ*  )+g(f)]<«  ♦  0(»j'1[||f"|Hlgll]) 

2nd  Step  Estimation  of  hj  and  h2 . 

From  (2.6).  (2.8)  and  (2.14), 

hj  -  J*  K^(f)tiij2f(f)+g(f)]<lf  *  k2f  (1 ) . 

where 

Kj  It )  *  -jfcosh  rj  (l  -(  J  t-cos  v(l-()]  -  ■|kji7_1(slnh  7(l-f)-sin  17  ( 1  -<■  ) ) 
Integration  hy  'arts  twice  for  f  yields 

Klrj  (Oiv2r(()<l(  -  f*  K^U)if"U)df  -  kf f ( 1 ) . 

where 

Klv((  )  -  \  (cosh  17  <  1-f  )-cos  t|(l-())  -  -j  kjT)-1  [slnh  ij(  1-f  )  »sin  v  (  W  ) ) 
We  get 

KinU>  -  £»f1e,,(i)-ik2)e'’,<  *  0(1). 

Klv(()  ’  jif1e’,<*»-lk2)e'’,t  *  0(1). 

Therefore 

h.  *  VM  (ii- Ik2  if1  e^f  [U"U)»gU)]d(  ♦  0<||HMIell )  ■  (215) 
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SiaiJarly,  from  (2.6)  and  (2.8). 

h2  =  Jg  K2Tj(f)[in2f(n+g(f)]df  -  k|f *  ( 1 ) 

itrt 

*2,(0  *  1  [sinh  v  (  W  )+sin  ij  ( l-f  )  ]  +  -|k|[cosh  77  ( l-(  )  -cos  ij(l-f  ) ) . 

topeating  this  same  Integration  by  parts  procedure  twice  more  for  f, 

m  get 

Ip  ^(Oi’J^fOdf  -  K2lI(f)if"(Odf  - 

to 

K2ij(f)  “  \v~1ev(l  +  iklv)e~v(  *  0(1) 

K2tj(0  -  J7i'1e’j(l*lk|n)e'’,f  ♦  0(1). 

N  get 

h2  *  ^■le’,U  +  lk|’»)J^  (lf"«  )-g«  ))d(  ♦  0(||  f"IM|  g||  ).  (2.16) 

3rd  Step  Estimation  of  Aj.Aj.A-j  and  A4. 


We  first  write 


Chen,  et  al 

11  i  i 

ii  -i  -i 

M  ‘  Je7!  (-i»j-lkj)e^  (-rlk^)e^  (l7j-ik*)e'1,J 

<1  +  i^k|)e^  (-l-ijk|)e^  ( l-l^k|)e^  (-i^k|)e‘l7l 

So 


Further,  write 

'•  *  "ll  "18 

Mj1  =  (det  Mj)-1  *  *  "21  "22 

'  *  "31  "32 

*  '  "41  "42.  . 

From  the  evaluation  of  cofactors, 

"ll  -  (l+l)[(l-»ik|)e~I,J  ♦  l(l-M|k|)8ll?  +  (l+l)(l-lijk|)e"’,)f 
~>il2  *  (l-l)[(’»-kf)e‘1’'  -  1  (ij»kj  )el7l  ♦  ( l-l )  (T/  +  ikj  ] . 

Let 

Hv)  •  (Tj-lkj)Mjj  +  (l  +  ii7k|)p12 

-  ( 1  +  i  )■>){  [••2k|’i+(  1*  1 )  ( l+kjk|)-21kjT7’J 
*(21k|»f«-(l*-l)(l+kjk|)'*-2kj»j*1  Jelll) 

♦  0(v2e~v) 

The  term  In  braces  above  satisfies 

l<  M  >  |  Bracket  2|  -  (Bracket  l| 

— *  2(l*k?k|),  as  -r)  ■»  ~ 
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if  c'v((ir"(t)>eu>i'Ui  *  et*  ^[iini*iiRin>. 

J  0 

Thus 

j  1 

Al’,2,K1I0e^!lf"(t),8(t)ldi  ♦  O^e^llir'IMIgll]).  (2.20) 

As  for  A2>  we  have 


A2V2  ”  Dl1  ^21^1*^22^2^ 

(2.21) 

where 

-m21 

=■  ( 1--1 )  ( l^-XTiklje17  ♦  2(l-77l<|)e''iT,  *  ( 1  +  1 )  ( l-lijk| )e'n 

(2  22) 

M 

22  *  ( 1-1 )  (»7  -  ikj  )en  -  2i(7j-k2)e'1’’  -  (  l*i )  (ij-ikf  le-’’ 

(2.23) 

By  (2.15), 

(2.16),  (2.22)  and  (2.23),  the  dominant  terms  In  /. i21 

hj  and 

P22h2  are 

that  they 

0(ije2v).  But  their  coefficients  In  li2ihi  *  /i22h2 
cancel  out.  We  get 

are  such 

Az*»2  -  O'1  •  0(nen[||f"lhll8ll  ] )  -  0(||r||*||g||) 

(2.24) 

Similarly, 

we  can  show  that 

a3t,2  -  0(|| f" INI ell) 

(2.25) 

a4t,2  -  O(iirilMlell) 

(2.26) 

Pinal  Step  Estimation  of  II  w** II  »||  vil  . 


By  (2.14)  and  (2.20),  we  have 


7T  t 
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w"  ( x )  --  w“  |x)  ►  w“  ( x ) 

*  (-  ^'1eT’x|g  e^[lf"(n*e(n)df  ♦  0(^~ 1  [|| f "U  +|le|! ] > } 

♦  Aj>;2(;’,x  +  A2(  tj2)o1,,x  ♦  A3(  '7,)2e  ’'x  -  •  |->;2)e  ,,x 

1 

-  {0(n_1I||f"|l-llell ] )  ♦  e7”1  •  0(»75e  T’[||f"IMIg||])} 

*  (A^e1”*  *  A3,2e  ’»x  -  A,,2.-1’*). 

In  the  first  parenthesized  term,  the  I2-norm  of  evx  is  of  order  of 
Magnitude 

[J*  (e’,x)2dx]1/2  «  [i^(e2n-l ) 1/2  *  0(7,‘1/2en)  . 

hence 

[f  |first  parenthesized  term|2dx]1/2  «  0(||f"|l*llgll )  • 

J  0 

The  second  parenthesized  term  also  satisfies  the  above  bound,  because  of 
(2.24) ,  (2.25) ,  and  (2.26)  . 

Hence 

||w»||  <  C(liriMkll)  -  0(||f"|MI*ll).  for  7,  large.  (2.27) 
For  v,  by  (2.4)  we  have 

IM|  <  |X|||w||  ♦  ||f ||.  (2.28) 

We  want  to  show  that 

IM2IMI2  S  c ( || w" II 2 *11  f" II 2 HI fll 2 *-11  ell 2 ) .  (2.29) 

for  some  constant  C  >  0  Independent  of  X 

Consider  (2.2.1),  with  o>4  -  1  and  x  «  in2,  and  use  k2,  k|  for 

2  and  k| . 

We  get 


Multiply  (2.2.1)  by  w(x)  and  integrate  by  parts  twice. 
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*  -<Xf-g,w> 

-  -<  f , Xw>- <g , w> . 


W»*(1JW(1)  V(I)K'|1)  -  II W" II 2  -  V’lMI2 


From  (2.2.3)  and  (2.2.4).  we  get 

i7,2k2|  w(  l )  |  2-k2f  ( 1  )w(  1 )  -b)k2|  w*  ( 1 )  |  2-k2f'  ( 1  )w*  ( 1 )  -||  w" ||  2-i,4||  w||  2 


«  -<f , Xw>~<g, W> . 


Hence 


>;4||w||2  >  Re{||w”||2+<f ,Xw>-<g, w>-k2f ( l)w(l)-k|f'(l)w'(l)} 

S  ||w"||2  -  -^)Jw|)2+2||f||2-  |j|gr)j2  -  |||„|)2  -  C[||f"||2-||w"||2]. 


where  we  have  applied  the  Poincare  inequality  and  the  trace  theorem: 


I  f  ( l )  I  2  -  |f  ( i )  I  2  <  c||f"||2 

|  w(  1 )  I  2  -  |  w*  ( 1 )  I  2  <  C||  w*’ || 2  . 

Therefore  (2.29)  follows  for  n  sufficiently  large. 

By  (2  27)  and  (2.29),  we  have 

||v||  <  C ( || f " || * || f || * || w" || ♦  || g|| )  (2.30) 

S  C (Ilf II *11  Kill 


as 


Ilf II  <  C||f|| 

and  (2.27)  holds. 

Combining  (2.27)  and  (2.30).  we  have  proved  (2.5)  for  |X| 
sufficiently  large,  X  »  1  u,  u  6  R.  So  Lemma  3  lias  been  proved.  Q 

THEORF.M  4  Let  k2  i  0  and  k^  >  0  in  (1.1)  Then  the  uniform 
exponential  decay  of  energy  (1.2)  holds. 
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Proof :  In  order  to  apply  Theorem  1.  we  need  to  verify  that  assumptions 

(1.6),  (1.7)  and  (1.6)  arc  satisfied. 

He  note  that  (1.6)  is  satisfied,  because  A  is  dissipative  and 

||exp( tA)||  <  1. 

The  verification  of  (1.7)  and  (1.8)  is  accomplished  if  we  can 
verify  merely  (1.7): 

(AI-A)-1  exists  for  all  A  «  io.  o  e  R.  (2.31) 

because  by  (2.31)  and  Lemma  3. 

KH  +  Hvxl!  -  C*  (||f"IMIg|l).  VA  -  io.  o  6  R. 


where 


C'  -  max(C,C").  C  as  in  (2.27)  and  (2.30) 


C" 


max  ||  (A  I -A  )-1||  ,  for  some  B,  sufficiently  large, 
|X|SB,  2 


A  ■  io,  o  €  R, 


To  show  (2.31),  we  assume  the  contrary  that  o(A)  f)  (i“l<v  6  R)  *  *■ 
By  Lemma  2.  o(A)  consists  solely  of  isolated  nonzero  eigenvalues 
Without  loss  of  generality,  let 

XQ  6  o(A).  AQ  -  It,2.  Vq  €  R.  t,0  *  0. 

Then 


has  a  nontrivial  solution  (Wq.Vq)  6  0(A).  Explicitly.  ( wy . vQ ) 
satisfies 
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i^o  -  v0  *  0  on  fo.l] 

f  4  j  2 

W0  *  iv0v0  “  0  0,1  t0'1) 

wQ(0)  *  0 

Wq (0)  >  0 

WqU)  -  kfv0(  1 )  -  o 

w”(l)  *  k|vp ( 1 )  =■  0 


Letting 

inn  t 

w(x,t)  -  e  Wq ( x ) , 
we  easily  check  that  w  satisfies 


Also,  the  energy 


“tt  *  "xxxx  -  °- 


/  nwxx(x.t)|2  ♦  |wt(x.t)|2]dx 


is  constant,  thus 


d 

ar 


f  ll»KX(x.t)l2  ♦  |wt(x.t)|2]dx 


*  0 

-  2Re(wxx(x.t)5xl(x.t)-wxxx(x.t)5t(x.t)]|^> 

-  -2(k2|wt(l,t)|2-k||wxt(l,t)|2) 

Because  k2  >  0  and  k2  >  0.  we  deduce 


|wxt(l  ,t)|  -  ’>olwoml  *  °- 
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•|wxxx<l.t)|  -  |Wg(l)|  -  0.  If  k*  .  0: 

|»xxx(l.t)|  -  l-Sdll  -  0.  |wt(l.l)|  -  ’*0 1  w0  ( i ,  I  *  °-  lf  kl  >  0 
Thus  Wq(x)  is  a  solution  to  the  boundary  value  problem 


Wg4)  -  Vq*0  -  0  on  [0.1] 
w0(0)  -  0 

"a (0)  -  0 

0  (2.32) 

*Q ( 1 )  *  0 

Wg(l)  -  0 

w^'(  t )  -  0 

Write 

w0(x)  -  A0Jcos  h0(xQ-l  l+Aggsin  v0(x  l)*A03cosh  »)0(*  1 )  ►A^sinh  nQ(x  l) 
Then  the  five  boundary  conditions  in  (2.32)  require  that 


cos  V0 

-sin  t,0 

cosh  Vq 

-slnh  v0 

A01 

sin  -n0 

cos  i)0 

-slnh  v0 

cosh  n0 

A02 

0 

1 

0 

1 

A03 

-1 

0 

1 

0 

A04 

0 

-1 

0 

1 

has  a  nontrivial  solution  (A01 . Aoz. A03> AQ4 ) .  However.  It  Is  easy  to 
check  that  the  matrix  in  (2.33)  has  rank  4  for  any  v0  e  K.  vQ  *  0.  a 
contradiction . 

Therefore  the  proof  of  Theorem  4  is  complete.  Q 


3 .  ASYMPTOTIC  ESTIMATION  OK  E 1  OKNKKKQOKNCl ES 

From  the  graphs  in  [2],  one  notices  that  at  low  frequencies 

eigenvalues  of  the  damped  operator  A  seem  to  exhibit  a 

"structural  damping"  [3]  pattern.  Does  the  structural  damping  pattern 
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continue  into  high  frequencies ,  or  is  it  only  a  low  frequency  phenomenon, 
for  beans  with  boundary  dissipation?  To  answer  this  one  must  examine  the 
asymptotics  of  etgeufrcquencles . 

The  work  of  asymptotic  analysis  was  first  done  by  P.  Rldeau  in  his 
thesis  (8]  (cf.  the  acknowledgement  at  the  end  of  the  paper)  Unaware 
of  his  results,  we  had  carried  out  the  analysis  independently.  We  feel 
that  it  is  of  significant  interest  to  include  the  work  here  as  it  will 
make  the  study  in  this  paper  more  complete,  and  only  a  minor  effort  is 
required. 

Let  (0(x) . ifr ( x ) )  be  an  eigenfunction  of  A  belonging  to  the 
eigenvalue  X(*0).  Then  by  (2.2),  setting  f(x)  »  g(x)  •  0  and  wx  *  0, 

we  see  that  0  satisfies 


«V4,(x)  ♦  X20(X)  -  0 
0(0)  -  0‘ (0)  -  0 

0"(i)  -  Xk2*(l)  -  0  (3.1) 

0" ( 1 )  ♦  Xk2*' (1)  -  o 


To  simplify  notations,  we  consider  the  following  eigenvalue  problem 


d(4,(x)  ♦  X2*(x)  •  0 

0(0)  -  0 
0'  (0)  •  0 

0-'(\)  -  Xk2*(l)  •  0 
•" ( 1 )  ♦  Xk|#' (1)  »  0 


Noting  that  the  following  correspondence 


A 

«2kj  in  (3.1) 
«2k. 


x 

k2  in  (3  2) 


(3.3) 


is  in  effect. 

The  boundary  value  problem  (3.1)  has  a  nontrivial  solution  if  and 
only  If 


Euler-Bernoulll  Bean  Equation 


(3.4) 


1 

1 

1 

1 

M^X 

pVx 

pVx 

3 

3 

3 

3 

(,«3x2  k2X) 

(,i9X2-Xk2) 

(/j,5X5-k2X) 

(,i2,X2-k2X) 

.  e^x 

.  e^X 

. 

• 

3 

3 

3 

3 

(A.2X*k2/iX^) 

(/A<p3k2X^) 

,  10,  5.,  2.  2. 

(/j  kjX  ) 

(p]4xVk2X^) 

.  e/V* 

.  e>*Vx 

.  e*1  A 

•  e^A 

-  0. 

where  /i  is  the  eighth  root  of  unity,  exp(ie/4).  The  derivation  of  the 
above  is  Identical  to  ( 2 . 10 ) ••  ( 2 . 12)  . 

Evaluating  this  determinant  yields  the  transcendental  equation 

2/2k|x{  I.-VW -eV2X  ,ie  V2X 

V^8(l  kfk| ) *2 ( 1  ■►k2k|)c'v'2X*2(l*k2k2)c  lv/2X  (3.5) 

♦2  ( 1  *k2k| )  </  ^  -2  < 1  *k2k2 )  e  ^2X  j 

^lk2{-ie-^-eV»McV^V2X}  ,  Q 

Now,  write 

X  -  |X|ei0.  (3.6! 


As  the  closed  right  half  plane  does  not  contain  any  eigenvalues,  and 
because  In  (3.5).  x  is  symmetric  with  respect  to  the  real  axis,  we  need 

only  consider  5j  <  9  <  n  In  (3  6).  We  will  actually  first  consider 

j  <  0  <  n  -  4 .  for  any  i  >  0  sufficiently  small  (3  7) 

The  case  of  0  -*  e  will  be  considered  in  (3.111-13.12) 

S  i  nee 


/X  *  | X| 1/2exp(  10/2)  -  | X | 1/Z[cos(0/2)*i  sin(0.-2)|. 
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we  see  that 


>V2X  =  e-/2pTcos(e/2)e-i/2]T[  sin(e/2) 


.V2X  .  -/2TX]sln(e/2)  i/2Txlcos(e/2) 


are  0(c-r,/^")  for  some  y  >  0.  Thus,  from  (3.5) 

♦  e"lv'2X[2W'2k|x«-2(Uk^k|)v/X  -  l2/2k2) 

*  e^Xf  2v'2k|x*2y/X(  l+kjk|)+2/2kj  ]  *  0(  |  X|  exp(-7/TxT) ) , 
which  implies 

g/2[”x"r  (-sln(0/2)»cos(0/2))  _  _e- i/flTf  [cos  (0/2 )  »s  in(0/2  )  ] 


(  2  V2k?X  *2  ( 1  *k?k?  )/X -i  Z/lv}.. 

| - f - - - - ij  *  0  ( |  X  I  exp  ( -vrq 

2/2k|x  i-2(l  *kjk|  )yx  *^/2k 2 


)) 


But  (assuming  k?j  >  0)  the  term  in  braces  equals 


i  1  Hill  e-i0/2  1  0(|X|-1). 

✓2  k^IxT 

Thus  we  seek  X’s  satisfying 

gy2jTf (-siri(0/2)*cos(0/2)]  (3  01 

2  2 

«  [.j  -  U-i )  (1*klk2)  e- 10/21  p- i^Fpr)  (cos (0/2 )*sin ( 0/2 )  J 

I-  kg/JXJ  J 

♦  0  ( |  x  | ' 1 ) 

We  observe  immediately  that  0  — *  tt/2  as  |X|  — *  »,  since  the  LHS  of 
this  equation  would  decrease  to  zero  otherwise.  Furthermore,  the 
equality  can  be  satisfied  (up  to  higher  order  terms)  only  when  the  first 
term  on  the  KHS  is  a  positive  real  number.  Thus 
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e  -  tyjpTf  [cos  (0/2  )  *s  in  (0/2 )  1  ~  e-V2fXT  •  /2  ~  , 
2SJTJ  Z  ( 2n  -  |)e. 


or 


n  is  a  large  positive  integer. 


(3.0) 


2 

The  above  gap  of  0(n  )  for  eigenvalues  is  common  for  Euler-Bernoul li 
beams  with  energy  conserving  boundary  conditions.  Now  we  see  that 
Euler -Bernoul 1 1  beams  with  boundary  energy  dissipation  also  have  this 
property. 

One  checks  that  when  the  RHS  of  (3.8)  is  real,  its  modulus  is 


1 


(1 -kjk|) [coa(9/2)+sin(e/2)  ] 

y2k*/!TT 


*  0  ( I  X(  ■'  )  . 


This  in  turn  implies  that  the  exponent  on  the  MIS  of  (3.8)  must  satisfy 


✓Ikfj/fxT 


If  we  now  write  0  «  j  ♦  t.  e  >  0,  and  expand  to  lowest  order  in  c, 
we  have 


(l*k?k|) 

*  ~  — r-5-  ■ 

k||X| 

Now  suppose  X  »  (  *  Ijj  .  Then  9  *  tanMij/f)  and  |a|  *  ( f 2  *-tj2  ) 1  /2 . 
Expanding  tan’1  about  >i/(  *  ■»,  we  have 


_  L 

V 


or 


3 


•fc+cjarv  ,  ,v/ 


S’ 


A’.l 
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L  _  U*k2k2) 

”  =  k  fTP^'i572 


Hence 


( 


^(l^k2^).  as  |\|  -*  ». 
k2 


(3.10) 


By  (3.7),  the  only  remaining  case  to  be  considered  is  when  Q  ■*  n, 
l.e.  when  X  approaches  the  negative  real  axis.  We  write 

✓X  -  |X| 1/2(cos(0/2)*isin(e/2)) 

as  before,  but  this  time  was  assume  |n-e|  <  <5,  with  4  small  —  say 
0  S  4  <  n/8.  Then  one  easily  checks  that 


|ev2x|  <  c. 


|  e1«/2X)  $  c 

for  some  c  >  0  so  that  (3.5)  can  be  rewritten  as: 


(3.11) 


<2^2  lk2X+2y2(l+k2k|)-2v'2ik2}e'/2TxT(sln!0/2)'lc°s(0/2) ) 

-  (2v'2k2X*2v'X(l+k2k2)*2v'2kj)el/2l>ll  (cos(9/2)*lsin(©/2) ) 

♦  O(X) .  (3.12) 

However,  for  9  in  the  range  of  interest.  sin(0/2)  >  2cos(0/2)  and  in 
particular,  sin(0/2)  >  0.5.  Thus,  the  modulus  of  the  L.H.S.  of  (3.12) 
will  be  much  larger  than  that  of  the  R.H.S.  (for  |X|  large)  so  this 
equation  has  no  solutions  if  |X|  is  large 


THEOREM  5  Let  X  =  (  *  in  be  an  e Igenf requency  of  vibration  of  the 
beam  equation  (1.1).  Then  for  |X|  large, 

r(2n-i)n,  2r_n  1/2 

|X|  ~  - g -  |S-I  ,  n's  are  large  positive  integers. 
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[»(EI)] 


1/2  ( 


{l+k2k|[m(EI  )]'*} 


(3.13) 


Proof:  Just  use  (2.3).  (3.3).  (3.9)  and  (3.10). 

By  (3.13),  the  eigenvalues  will  be  distributed  nearly  parallel  to 
the  imaginary  axis  at  high  frequencies.  Therefore  there  is  no"structural 
damping”  when  the  frequencies  are  high.  This  has  also  been  confirmed 
numerically  in  Rldeau's  thesis  [8]. 

We  note  that  when  k|  -  0  and  k2  >  0 ,  asymptotic  limits  can  be 
obtained  in  the  similar  way. 


4.  DESIGN  OF  PASSIVE  DAMPING  MECHANISMS 


The  following  is  a  (more  or  less  exhaustive)  list  of  combinations 
of  dissipative  boundary  conditions  for  an  Euler-Bernoulll  beam: 


-ElyxxU.t)  ■ 


-k‘yt(l.t) 


kj  >  0 


(4.1) 


EIyxxx* 1 ' 


-Eiyxxn.t>  -  *Syxt»>.t> 


x|  >  oj 


(4.2) 


yxd.t)  *  o 


-Eiyxxxd.t)  *  -kjyt(l  .t) 


kj  >  o 


(4.3) 


y(i.t)  -  o 

-Elyxx(l.t)  -  k22yxt(l.t) 


k|  >  o 


(4.4) 


EIyvvv(l.t)  -  -kjyt(l.t)  *  C1yv,(l.t 


1  *  xt  1 


k 2  >  0. 
k|  >  0 . 


-Ely^ll.t)  >  k|yxt(l.t)  *  c2yt(l.t) 
where  in  (4.5),  Cj  and  c2  are  real  constants  satisfying 


(4.5) 


(Cj-C2)a0  -  k2a2  -  k202  SO  V  a.p  6  R . 


(4.6) 


jgiy! .  1  .  j 


w- ■ 
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Note  that  the  boundary  conditions  (1.1.4)  and  (1.1.5)  correspond  to 
C1  ”  c2  “  0  *n  d-5)-  Obviously,  (4.6)  Is  satisfied  in  this  case. 

We  want  to  show  that  all  stabilization  schemes  (4. 1  )  —  (4.0)  can  be 

realized  in  practice,  at  least  by  designing  passive  dampers. 

As  (4.5)  seems  to  represent  the  most  complicated  case  among 
(4.1) -(4.5),  we  treat  it  here,  at  least  for  certain  special  values  of  Cj 

and  c2  (cf.  (4.7)  later).  The  other  cases  can  be  studied  similarly. 

The  following  damper  arrangement  gives  a  design  which  effects  the 

coupling  of  shear  (resp.  bending  moment)  with  velocity  and  angular 
velocity: 


a)  Inclined  Damper 


b)  Shortening  Velocity  c)  Damper  Forces  on  Beam 

of  Damper 

Shear(l,t)  =  -c,  v  sin  9 
d  s 

Moment(l,t)  «  c,  v  cos  9  h/2 
d  s 


Figure  2  Damper  arrangement  for  (4.5) 


Euier-Bernoulli  Beam  Equation 


A  single  damper  (cf.  Figure  2a)  is  attached  to  the  lower  end  of  the  beam 
at  an  inclination  angle  0  with  respect  to  the  horizontal.  Using  the 
velocity  at  the  end  of  the  damper,  vs,  and  the  associated  forces  shown 

in  Figure  2b  and  2c,  we  obtain 


Shear(l.t)  *  ~cdvs  s*n  ® 

Moment (l.t)  =  cdvs*cos  9)g  , 

where  cd  represents  the  damping  coefficient  associated  with  the  damper 
in  use.  As 


Js  -  yt(l.t)sin  0  -  yxt(l,t)|cos  9, 


we  get 


Shear ( 1 ,  t )  *  *EIyXxx(  1 ' =  *(cdsin2e)yt  <•  (|  •  edsin  9  cos  0)yxt 

2 

Moment ( 1 ,  t )  *  -E I yxx (l.t)  -  ( £-cdcos20)yxt  ♦  (-  l;'cdsin  0  cos  0)yr 


A  comparison  of  the  above  with  (4.5)  shows  that 


k2  =  cdsin20:  k|  =■  jj-cd  cos29 
Cj  «  -c2  *  gc^sin  0  cos  0, 


(4.7) 


thus 


(Cj-CgJa^  -  k2a2  -  k|fl2  *  -(kjO-kg^)2  <  0  V  a  ,fi  G  R 


so  (4.6)  Is  satisfied  and  the  boundary  conditions  (4.5)  are  dissipative. 
It  is  noted  that  when  0  «  n/2  (vertical  damper),  the  above  gain 
2  2 

constants  reduce  to  kj  »  cd  t  and  k|  *  Cj  =  c2  »  0,  cf.  Figure  3a. 

Similarly,  for  0  *  0  (horizontal  damper).  Hie  gain  constants  become 
2  2  2  • 

k g  *  cd  2h  /4  and  kj  *  Cj  »  c2  *  0,  as  shown  in  Figure  3b. 
Consequently ,  the  boundary  conditions  ( 1 .  1 .4 ) - ( 1 .  1 . 5)  can  be  realized  as 
in  Figure  3c,  with  k?  -  cd  .  and  k2  =  c.  2h2/4 


Chen,  et  al 


I  ca  yt(1*E) 

Shear(l.t)  »  -cd  y^d.t) 
a)  Vertical  Damper 


cd  yxtn.t)h/2 


Moment(l,t)  =  y  (l,t)  h^/4 
b)  Horizontal  Damper 


x«  1 


c)  Shear ( 1 , t)  =-cdyc(l,t) 

Moment (1,C)  *  cdyxt(l,t)  ^  ^ 

Fieure  3  Damper  arrangement  for  (4.1),  (4.2)  amJ  ( 1 . 1 . 4 )  +  ( 1 . 1 . 5  ) 
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The  other  boundary  conditions  (4.1)-  (14)  can  be  realized  and 
designed,  respectively,  as  In  Figures  3a.  3b.  4  and  3 

The  method  of  estimation  which  we  have  developed  in  this  paper  and 
Huang's  theorem  (Thm.  1)  can  be  applied  to  study  exponential  stability 
for  all  of  these  boundary  stabi  1  iz.at i on  schemes 
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